Flow equations for the one-dimensional Kondo lattice model: Static and dynamic 

ground state properties 
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The one-dimensional Kondo lattice model is investigated by means of Wegner's flow equation 
method. The renormalization procedure leads to an effective Hamiltonian which describes a free 
one-dimensional electron gas and a Heisenberg chain. The localised spins of the effective model are 
coupled by the well-known RKKY interaction. They are treated within a Schwinger boson mean 
field theory which permits the calculation of static and dynamic correlation functions. In the regime 
of small interaction strength static expectation values agree well with the expected Luttinger liquid 
behaviour. The parameter K p of the Luttinger liquid theory is estimated and compared to recent 
results from density matrix renormalization group studies. 



I. INTRODUCTION 



The fascinating subject of heavy fcrmion physics in 
rare-earth and actinide systems has been a challenge for 
theoretical and experimental investigations for decades 
The intriguing properties of these materials are far 
from being understood and still give us a lot of puzzles 
to solve. Theoretical studies of the heavy fermion ma- 
terials are based on several models like the periodic An- 
derson model (PAM) 0. Another generic model is the 
Kondo lattice model (KLM) which describes a noninter- 
acting electron gas coupled to localised spin moments via 
a Heisenberg spin interaction. The Hamiltonian reads 



Sk c!c k0 



r af3 Ci/3, 



(1) 



ka 



i af3 



where Sk — — Ylij Uj e lk ^ Ri ^ R ^ is the dispersion relation 
for the electrons on the lattice, being the hopping 
integrals. The parameter J is the exchange integral of 
the local spin interaction, the so called Kondo exchange. 

We want to consider here the one-dimensional case 
which has been the subject of numerous numerical and 
analytical investigations. Numerical studies were based 
on the Quantum Monte Carlo (QMC) method [|[, exact 
diagonalization (ED) studies [j,Jal, the density matrix 

renormalization group (dmrg) [aaSiO 

or the nu- 
merical renormalization group (NRG) method [n|. Ana- 
lytical approaches comprised the bosonisation technique 
[H, [H[ or the renormalization group (RG) theory [l4 |. 

The phase diagram of the one-dimensional KLM as a 
function of the Kondo coupling J and the band filling 
n c of the conduction electrons is quite accurately known. 
In higher dimensions the KLM is believed to show the 
well-known Doniach phase diagram [151 ]. In contrast to 
the latter the one-dimensional model does not exhibit a 
magnetically ordered phase in the parameter regime of 
small interaction strengths J. In this parameter regime 
the corresponding phase diagram is governed by a para- 
magnetic metallic phase |4|. There, the model is assumed 
to belong to the universality class of the so called Lut- 
tinger liquids [l6[ which possess gapless charge and spin 
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FIG. 1: Phase diagram after [12 



excitations resulting in an algebraic decay of the corre- 
sponding correlation functions. The asymptotic form for 
density-density- and spin-spin-correlations are [17| 



(Sn(x)Sn(O)) 



(S(x)-S(O)) 



K„ 



(tt x) 2 



Ax cos(2k F x) x^ 1 ^ Kp 
A 2 cos(4:k F x) x~ 4Kp 



Bi cos(2fcj?a;) x~ 



-1-K a 



(2) 



(3) 



The parameter K p is a model dependent constant which 
determines the low-energy physics. Apart from the para- 
magnetic metallic phase for small J ft the phase diagram 
further comprises a ferromagnetic ordered phase for large 
J and a spin liquid insulator phase at half-filling, n c = 1. 
There are two limiting cases in which the ground state 
has been proven to be ferromagnetic Q ■ Firstly, the limit 
of vanishing electron density n c — > 0, secondly the case of 
infinite coupling strength J ft — » oo. The situation at half 
filling is special in the sense that it exhibits finite gaps 
for spin and charge excitations at any finite coupling J. 

The KLM can be understood as an effective Hamil- 
tonian of the above mentioned PAM. It is connected to 



2 



the PAM by a Schrieffer- Wolff transformation [18]. This 
property naturally raises the question whether the lo- 
calised spins in the KLM participate in the formation of 
the Fermi surface, or in other words: Does kp = n c 7r/2 
or kp = (n c + l)7r/2 hold? The size of the Fermi surface 
can be read from the positions of singularities in certain 
correlation functions. Recent results seemed to confirm 
the picture of a small Fermi surface with kp = n c n/2 [7|. 
However, a more careful analysis which has recently been 
performed by Shibata et al. supports a large Fermi 
surface. 

In this paper we shall apply the analytical method of 
continuous unitar y tr ansformations (flow equations) pro- 
posed by Wegner [20| and Glazek/ Wilson [2l( to the one- 
dimensional KLM. It was first app lied to this model in 
arbitrary dimensions by Stein [221 ]. He derived an ana- 
lytical expression for the RKKY interaction. 

In Sec. |TT] we shall give a short introduction into the 
flow equation method. In Sec. IHII the method will be ap- 
plied to the one-dimensional KLM. By integrating out the 
Kondo coupling between the conduction electrons and lo- 
calised spins we arrive at a decoupled system of a renor- 
malized noninteracting one-dimensional electron gas and 
a renormalized spin chain. In the latter the spins inter- 
act via an effective spin exchange. Within the frame- 
work of the flow equation method it is straightforward to 
find expectation values and correlation functions, if the 
eigenvalue problem of the effective model is known. In 
Sec. IIVI we shall show how the method can be used in 
order to verify the expected characteristic behaviour of 
a Luttinger liquid. Previous investigations of the one- 
dimensional KLM have mainly focused on static proper- 
ties like the momentum distribution or spin and charge 
correlation functions. In this work we shall put special 
emphasis on the investigation of dynamic properties and 
extend already existing results for the dynamics. 



II. FLOW EQUATION METHOD 



equations by using the expression 

dU(l) 



'7(0 = 



dl 



U\l) 



(5) 



for the antihermitean generator = — 77^ (Z) of the uni- 
tary transformation. The differential equation for the 
Hamiltonian takes the simple form 



dH{l) 
dl 



= W),H(l)} 



(0) 



The generator has to be suitably chosen. Wegner's ap- 
proach starts from a decomposition of the Hamiltonian 
into an unperturbed part Tio, whose eigenvalue problem 
is assumed to be known, and a perturbation H±. Weg- 
ner's generator is given by 



r,(l) = [H (l),W)} 



(7) 



which is simply the commutator between the unper- 
turbed part Ho(l) and the perturbation 7ii(l). This gen- 
erator integrates out all interaction terms except for pos- 
sible degenerations [2(J. It finally leads to a diagonal or 
block-diagonal effective Hamiltonian. 



III. FLOW EQUATIONS FOR THE KONDO 
LATTICE MODEL 



We now turn to the derivation of the flow equations for 
the parameters of the Hamiltonian. With this in mind 
we proceed as follows. Firstly, we give the flow invariant 
Hamiltonian which includes new generated, effective in- 
teractions. The flow invariant Hamiltonian then leads us 
to the specification of the generator. Thereby we shall 
introduce some of the necessary approximations within 
our approach. 



To begin with we would like to sketch the concept of 
the flow equation method which was independently de- 
veloped by Wegner [20] and Glazek/Wilson [U in 1994. 
Since then the method has successfully been applied to 
a great number of problems, e.g. the electron-phonon- 
problem (23j . one-dimensional interacting fermion sys- 
tems [lij or the spin-boson-problem (25j . 

The basic idea of the flow equation method is the ap- 
plication of a continuous set of unitary transformations 
to a given Hamiltonian 



H(i)=u(i)nu i (i). 



(4) 



Here I means the continuous flow parameter. The pur- 
pose of this procedure is that one wishes to diagonalizc 
or at least simplify the Hamiltonian. Thereby the pa- 
rameters of the Hamiltonian become renormalized. This 
treatment is translated into the language of differential 



A. Flow equations for the Hamiltonian 

The first step in deriving the flow equations is the de- 
termination of the generator r/(l). In a first step we wish 
to integrate out the Kondo coupling between the conduc- 
tion electrons and the localised spins, so the most simple 
generator is 



ikq a/3 



■<r a p c{ a c qP e l( - k ~ 



q)Ri J 



(8) 



where the coefficients rj^ q (l) are still unspecified. They 
depend on the concrete choice of the generator. Weg- 
ner's approach starts out from the generalised form of 
Eq. ([7]). If we take only the conduction electrons to be 
Ho, we obtain rjk q (l) = (£& — £ q )Jkq(l) for the coefficients 
of the generator. The commutator between the generator 
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([8]) and the Hamiltonian JT]) gives rise to new, effective 
interactions. Using Wegner's approach they enter the 
generator and are eventually integrated out. We shall 
introduce a more general form of rj(l) below. 

I 



where : X : denote operators resulting from a decoupling 
scheme which we shall discuss later. 

Before we proceed let us take a closer look at equation 
©. The first line represents the block diagonal part of 
the Hamiltonian since electron and spin operators are 
decoupled. It contains a complicated RKKY-like spin 
interaction term between the local moments. The second 
and third line comprise the nondiagonal or interaction 
part. Aside from the Kondo coupling we get interactions 
between the local moments which are either symmetric 
or antisymmetric with respect to interchange of the sites. 
Correspondingly, the first one couples to the electronic 
charge density, whereas the second one couples to the 
electronic spin density. We restrict ourselves to these 
terms because they are the most important ones in the 
regime of small interaction strength J. That way the 
above Hamiltonian becomes flow invariant and Eq. ([9]) is 
valid for all flow parameters I. For / = it represents the 
initial Hamiltonian (fT]). This implies the following initial 
values of the parameters 

£fc(Z = 0)=£ fe , J kq {l = 0) = J 

I 



and the prefactors T] kq (l), tj^ (I) and t]^ (I) are deter- 
mined by Eq. Q. 

After the transformation, i.e. in the limit I — > oo, 



In order to see what kind of effective interactions 
emerge, let us commute the initial Hamiltonian (fTJ) and 
the generator of Eq. ([8]). After some calculation we ob- 
tain the following Hamiltonian 



(9) 

I 

XkqQ = 0) = 0, M kpq (l = 0) = 

D kpq (l = 0) = 0, E c (l = 0)=J2 £k71k - ( 10 ) 

k 

The prefactor of any operator term of Eq. © con- 
trols the strength of the respective operator. Within the 
framework of the flow equation method they are deter- 
mined by corresponding differential equations. With the 
choice of the generator rj(l) we can control which of these 
terms are kept and which are to be vanished. Since the 
aim of our renormalization procedure is a blockdiagonal 
Hamiltonian in which electron and spin operators are de- 
coupled, we have to remove all terms describing interac- 
tions between electron and spin operators. The generator 
r)(l) of the continuous unitary transformation has to be 
chosen appropriately. 

With this in mind we can now write down the gen- 
erator r)(l). Using Wegner's approach we have to take 
into account the generated, effective interactions. The 
generator reads 



11) 

I 

only the first line of Eq. ([9]) remains. It represents the 
diagonal part 7io(l). This effective Hamiltonian can be 
used to easily calculate physical properties. The nondi- 



11(1) = ^2zk{l) ■ c\ a c ka : +^ ^Xk q {l) : S k - q ■ S 



kq 



q—k 



-E c (l) 



kq af3 



kpqa 



-2 ^2 iD kpq (l) (S fe _p x S p _ g ) ■ (T a p : c\ a c qj3 : 



AN 2 



kpqafl 



H e (l) + H s (l) + E c (l) + Hj(l) + Hm(1) + H D (l), 



(I) = rr 7 (l) + V M (l) + v D (l) = ±Yl Sk -i ■ ■ C L%3 : +7^2 E : Sk -P ■ V* ■ ■ c l 



kqa(3 



AN 2 



-ko^1° 



+ E Vkpqi 1 ) «(Sfe-p X S p _ ? ) • (T a 



kpqa 



(3 : C L%3 



kpqa (3 
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agonal part Ti.%(l) vanishes for / — > oo and the effective 
Hamiltonian TL := Ti.il — oo) then reads 



= 7i e + TLs + E c 



kq 



(12) 



In the following we shall denote all renormalized variables 
by a tilde. As Eq. (fl2|) tells us the effective model will 
consist of a one-dimensional noninteracting electron gas 
and a Heisenberg spin chain with renormalized parame- 
ters. 

We now have all ingredients needed to derive the flow 
equations for the parameters of the Hamiltonian. Before 
doing this we want to look at the approximations that 
have to be done. Firstly, we neglect interactions of or- 
der C(J 3 ) and higher in the Hamiltonian ©. Secondly, 
we decouple higher operator products in order to reduce 
them to those appearing in H(l). This gives rise to oper- 
ator expressions of the form : X :. They refer to fluctua- 
tion operators and mean either a normal order product of 
fermionic operators or a Hartree-Fock-decoupling scheme 
of spin operator products 



: c \a c ka '■ 



T C ko 



(13) 



' S 9 _fc I — Sk-q • S q -k — (Sk-q ' S 9 _fc). (14) 



The thermodynamic average will here be taken with re- 
spect to the effective model TL, Eq. (fl2|) . which de- 
scribes a decoupled system of a simple Fermi gas (elec- 
trons) and a Heisenberg spin chain with long-range in- 
teractions. These expectation values are therefore l- 
independent. The decoupling leads to a formal tem- 
perature dependence of the flow equations. Here we 
consider only the ground state properties, i.e. T = 0. 
For the sake of simplicity we introduce the abbreviation 
S(k — q) :— (Sfc_ g • S q ^k) for the spin correlation func- 
tion. One may also think of other expectation values 
like {Sk-q x S q -k) or (S q ). Since we consider the limit 
of small Kondo coupling J, the system is in the param- 
agnetic phase, where no symmetry is broken. Therefore 
these expectation values vanish. 

Evaluating the commutator between the generator (fTT|) 
and the Hamiltonian ([9]) we arrive at the flow equations 
for the parameters of the Hamiltonian. For the sake of 
clarity the /-dependence of all parameters is dropped. 
The electronic single particle energies £k obey the fol- 
lowing differential equation 



dl 



1 

2N 



(15) 



For the paramter \kq of the effective spin interaction 
we obtain the following flow equation 



dXkq , \ J T 

—jf - = {nk-n q )r) kq J qk . 



(16) 



The occupation numbers n k which enter the above equa- 
tion are again formed with respect to the effective model. 
The constant E c of TLq(1) follows 



dE c 

dl N 



lj^2( n k ~ n q ) S(k - q) f]( q J qk . 



(17) 



kq 



We restrict the renormalization of the effective interac- 
tion terms to contributions of order 0(J 2 ). Therefore 
both coupling parameters D kpq and Mj~ pq obey the same 
flow equation 



= \iVkpJpq + VqpJkp) - (s k - e q ) T$ m . (18) 



The first term is responsible for the generation of the 
effective coupling while the second contribution, which is 
always negative, ensures the vanishing at the end of the 
renormalization procedure. Finally for the flow equation 
of the Kondo coupling we find 

dJkq / \ J 

- -{£k-e q )r]i q 



dl 



^/^rP^^I ( 7 lkpJpq + VqpJpk) 
p \ / 

+ gJV Z^^kpDpkq + VqpDpqk) 

p 

+ y^d^fcqp Jpk + VqkpJpk) 

P 

~ iT/V ^ j JlkpDp,P+q— k,q + VqpDp^p+k-q.k) 



8N 

P 

-E 

8iV ^ 

p 



(Vk,p+q-k,pJpk + Vq,p+k-q,pJpk): (19) 



where we have taken into account correction terms up 
to order 0(J 3 ). Therefore we expect to find reasonable 
results only in the parameter regime of small coupling 
strength J/t. As this ratio increases further correction 
terms have to be included. The flow equations (fl~5| to 
(HHJ) represent a closed system of first order differential 
equations, whose solution can only be found by numerical 
integration. 



Here S(k — q) is the local moment's spin correlation func- 
tion which has to be evaluated with respect to the renor- 
malized Hamiltonian TL. It is therefore /-independent. 
As the effective model is not known before the end of 
the transformation we have to solve all flow equations 
self-consistently. 



B. Approximations for the effective model 

In the preceeding subsection we have derived flow 
equations for the parameters of the Hamiltonian. As to 
solve them we still need an analytical expression for the 
spin correlation function S{k — q). As it describes spin 
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correlations of the effective model, we are dealing here 
with a one-dimensional Heisenberg chain with long-range 
interactions whose exact solution is not known. Hence, 
we have to resort to further approximations. We stress 
here that this is the most crucial approximation within 
our approach because it strongly affects all renormalized 
quantities. Since the spin interaction is the result of the 
continuous unitary transformation it is not known until 
the transformation is completely performed. 
As our approach is only valid for small J/t, i.e. for the 
paramagnetic metallic phase with no broken symmetry, 
we use the Schwinger boson formalism to describe the 
spin system 26| . It preserves the rotational invariance of 
the spin Hamiltonian. The spin operators are expressed 
in terms of Schwinger bosons oti a and a\ a according to 



(20) 



where <7^ CT , stands for the Pauli spin matrix. Since the 
occupation number for bosons is not restricted, a local 



2S must be en- 



constrained of the form 
forced. 

We follow here the procedure of Trumper et al 
of Ceccatto et al. ual and introduce two fields 



A{j _ ^ (T QjtjCljfj 



= -A, 



and 



Bi. 



£4, 



B] 



or 



(21) 



(22) 



describing antiferro- and ferromagnetic correlations, re- 
spectively (a = —<t). This yields to the following Hamil- 
tonian 



■l.jXUf.Jhj 



(23) 



The expression M(0) stands for a normal order product 
of boson operators. The Hamiltonian is now biquadratic 
with respect to the Schwinger boson operators. We use 
a mean field theory in order to decouple the biquadratic 
terms. By using the mean field parameters (-By ) and 
(Aij) and replacing the local constrained by a global one 
we obtain a Hamiltonian which can easily be diagonalized 
via a Bogolubov transformation. Introducing new boson 
operators = u/. 4- iav^ a}_ ks we obtain 



Hs 



2^ 

qa 



L0 a 



(24) 



with ui(q) = \J (7b(q) — A) 2 — "j\{q) representing the en- 
ergies of the elementary excitations aJ CT of the spin sys- 
tem. Here the quantities ja (q) = § Y^Ri Jij(Aij) e lqRij 
and 7s(g) = \ Jij(Bij) e lqRij are used. The mean 



field parameters (-By) and (A^) and the Lagrange pa- 
rameter A have to be determined selfconsistently by solv- 
ing the corresponding saddle point equations. 

Finally we find an analytic expression for the spin cor- 
relation function which for T = reads 



S{q)T=o 



with Qi 



1 

AN 



■ tanh 



£(cosh[2(0 fc -0 fc+ ,)]-l) ) (25) 
fc 



Compared to methods like the Bethe ansatz for the 
nearest-neighbour Heisenberg chain the approximative 
Schwinger boson treatment discussed above has the ad- 
vantage that as many interaction terms as possible can be 
taken into account. With the approximation for the ef- 
fective model we are able to describe the one-dimensional 
KLM consistently within the framework of the flow equa- 
tion method. Any physical quantity we are interested in 
can be evaluated within the present approach. Especially, 
we emphasise that nothing has to be put in by hand. 



C. Expectation values and correlation functions 

We now turn to the calculation of expectation values 
and correlation functions. In this subsection we give the 
essentials for the derivation of certain important expec- 
tation values and correlation functions. We shall discuss 
the results in Sec. ITVl 

The retarded Green's function between operators A 
and B is in general defined as the following commutator 
or anticommutator relation 

Gab® - -id(t)((A(t);B}) = -iO(t){[A(t),B]±), (26) 

depending on the statistics under consideration. The 
thermodynamic average and the time-dependence have 
to be taken with respect to the full Hamiltonian. One can 
exploit the invariance of the trace under unitary trans- 
formations and obtains 



GAB(t) = -i6(t)((A(t);B)) n . 



(27) 



Now the thermodynamic average and the time- 
dependence are taken with respect to the effective model. 
According to the transformation of the Hamiltonian we 
also have to transform the operators. They obey a similar 
flow equation as the Hamiltonian 



dA(l) 
dl 



W),A{1)] 



(28) 



The commutation between r)(l) of Eq. and the 

electron operator Cka- leads to the following operator 
structure 



Cfco-(Z) = Q!Jfe(Z) Cka 



(29) 
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where we have taken only the correction terms into ac- 
count that couple to one local moment. The initial 
conditions of the parameters are ot k (l = 0) = 1 and 
lkq{l = 0) = 0. We transform the spin operator ac- 
cording to 



Sf (J) = + Cb,(0 vS'e^-^ clc qo 



kqcr 



(30) 



s? (i) = ms° + ^Y. M) ° s ? Sk ' q)Rl a ' c L<v 



kqa' 



q)Ri J 



kq 



(31) 



Here the initial parameters are (3(1 = 0) = 1 and C, kq (l = 
0) = 0. 

In order to derive the flow equations for the param- 
eters of the operator transformations we have to use an 
equivalent decoupling scheme as for the Hamiltonian. We 
finally arrive at the following differential equations 



Vkqlqk 



dlkq 
dl 



2 %k a k 



(32) 



(33) 



for the parameters of the electron operator transforma- 
tion and 



dp 
~dl 



J^^kqCkq n k {l-n q ) 



kq 



kq 



dl 



= -PVqk 



(34) 



(35) 



for the parameters of the spin operator transformations. 
We notice that the spin correlation function S(k — q) 
of the effective model enters the flow equation of a k (l) 
whereas the occupation numbers n k govern the flow equa- 
tion of 0(1). We restrict the flow equations for the correc- 
tion terms to first order contributions in the Kondo cou- 
pling. Going beyond this approximation could bring us 
up against the violation of certain summation rules which 
have to be fullfilled. We can combine the above equations 
to obtain flow invariant expressions. The expectation val- 
ues S(k — q) and n k are taken with respect to the effective 
Hamiltonian and are therefore ^-independent. We arrive 
at 



4(0 



l£s(fc-<zbi(0 = i 



and 



kq 



(36) 



(37) 



which displays the unitarity of the transformation. 

After determining the operator transformation we are 
now able to calculate static and dynamic correlation func- 
tions that characterise the ground state properties of the 
one-dimensional KLM. One of the most important quan- 
tities is the momentum distribution n(k) which reads 

n ( k ) = (Aa c k<y) = <5fc n k + — ll q S(k - q) n q . (38) 

For a Luttinger liquid we expect a continuous behaviour 
with respect to the momentum k and a power law sin- 
gularity at the Fermi momentum. The position of this 
singularity fixes the size of the Fermi surface. 

The static correlation function of the local moments 
Sff(q) indicates the phase transition from the paramag- 
netic phase into the ferromagnetic phase on increasing 
the Kondo coupling J. Within our approach it is given 
by 

S ff(l) = ( S 9 • S -q) 

= /3 2 S(q) + S(k-p + q) n k {\ - n p ). 



kp 



(39) 



We can also evaluate the static charge correlation func- 
tion C(q) and the static spin correlation function S cc (q) 
of the electrons. Their rather lengthy expressions are 
given in the appendix. 

The flow equation formalism allows us to calculate dy- 
namic quantities. The first quantities we look at are the 
one-particle spectral functions A±(k,Lu) of the conduc- 
tion electrons which measure occupied and empty states 
of the conduction electrons. 

/•OO 

A+{k,<j) = / dt(c ka (t)cl a )e iut 



a k (1 - n k ) 8(lo - e k ) 



\rl ^2 Ikq ( U p Vk+p-q - Vp Ufc+p-g) 2 



2N 2 



A_(k,u) 



x (1 - n q ) S(uj -s q -Q p - Lu k+p - q ) 

(40) 

dt(cl(t)c ka )e^ 



= a k n k S{uj - e k ) 

+ 2~^2 ( U p Vk+p-q - v P u k+p - q ) 2 

qp 

x n q 5(uj - e q + Lu p + w fe+p _ 9 ). 

(41) 

Here u k and v k are the coefficients of the Bogolubov 
transformation used to diagonalize the Schwinger bo- 
son mean field Hamiltonian (JMJ). The spectral func- 
tions A±(k,ui) comprise two contributions. The first 
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term (~ a 2 ) represents a coherent quasiparticle excita- 
tion. The second term is an incoherent background. It is 
important to note that the elementary excitations of the 
spin system of the effective Hamiltonian Cu q enter the lat- 
ter contribution. The electronic density of states defined 
by 



with G(k,Lu) being the electronic Green's function, can 
also be calculated. 



S k 



(42) 



Another important quantity is the dynamic spin struc- 
ture factor Sff(q,Lo) of the local moments 



/DO 1 
*(S g (t).S_ g )e*-* = — ^/3 2 



'p+q - V p U p+q ) 2 8{UJ -Ulp — UJp+q) 



2 ^ 

Ckp( u p' v p'+k-p+q - W+fc-p+?) 2 n k (l-n p ) 



kpp' 



(43) 



The first line describes only the spin excitations of the 
effective model in terms of Schwinger bosons. The sec- 
ond line of Eq. (|43]) results from the coupling of the local 
moments to electronic particle-hole excitations of the ef- 
fective Hamiltonian. In addition, we can also calculate 
the dynamic spin structur factor of the conduction elec- 
trons S cc (q,uj) which is found in the appendix. 



IV. RESULTS 

After having derived theoretical expressions for various 
correlation functions from the flow equation approach we 
now turn to present the outcome of the numerical solu- 
tion of the flow equations US]) - |Q]5J) and (02) - © ■ We 
start with the result for the parameters of the Hamilto- 
nian and subsequently show our findings for static and 
dynamic correlation functions. We shall show to what 
extend the statics reflects the expected Luttinger liquid 
behaviour. We also clarify the possibility of the approach 
to describe the quantum phase transition on increasing 
coupling strength. 



Parameters of the Hamiltonian 



In order to solve Eqs. (| 1 5|) - (fig)) we used a Runge 
Kutta algorithm. The complexity of the differential equa- 
tion restricted our system size to N — 120. Remember 
that the spin correlation function S(k — q) which enters 
the flow equations has to be calculated with respect to 
the effective model (fT2"j) . Therefore the parameters of the 
Hamiltonian had to be determined self-consistcntly. 

The spin correlation function S(q) of the effective 
model plays an important role. We therefore start our 
discussion with S(q) which is shown in Fig. [5] The main 
feature is the dominant peak that shows up exactly at 



2 
c55 




FIG. 2: Spin correlation function S(q) of the effective model 
for two different fillings n c . 



the wave vector q = 2kp = n c n, where hp is the Fermi 
momentum of the conduction electrons. As we shall see 
later the pronounced structure has severe consequences 
for other quantities that are related to S(q). The pro- 
nounced peak is due to the special excitation spectrum 
of the Schwinger bosons. The other main property of 
S(q) is the vanishing ferromagnetic component (q = 0) 
which can easily be understood from Eq. (f2"5)) . 

The elementary excitations u) q of the spin system of 
TC are shown in Fig. [3] They exhibit a small but fi- 
nite gap at q = k c F mod ir. This small gap is respon- 
sible for the strong peak in S(q). It manifests the ro- 
tational invariance of the ground states and is an arti- 
fact of the Schwinger boson approach as we are dealing 
here with half integer spins (S — 1/2) which may have 
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FIG. 3: Elementary excitations u) q at n c — 2/3 and J/t = 1.5. 



a gapless excitation spectrum. Nevertheless, within the 
Schwinger boson approach a vanishing gap would give 
rise to a ground state with broken symmetry that con- 
tradicts the assumption of a rotational invariant para- 
magnetic phase. However, the important point is the 
position of this gap. It determines the maximum of the 
spin correlation function which turns out to be at the ex- 
pected position. Therefore wc assume that a description 
in terms of spinons should not change these results de- 
cisively. Looking at Eq. (|2"5")) we see that always pairs 
of excitations enter the equation for S(q) so that the 
maximum of the spin correlation function is found at 
q = 2kp — n c n. 

At this point we add that we found solutions for the 
saddle point equations of the SBMFT only in the param- 
eter regime 1/2 < n c < 1. The case of half filling is 
special in the sense that there exists a gapped spin liquid 
phase. It remains an open question whether the present 
approach can also be used to describe this phase. Below 
n c = 1/2 the dominance of the ferromagnetic compo- 
nents in the RKKY coupling Jij prevents a solution of 
the saddle point equations of the SBMFT. 

Finally we discuss the renormalized electronic single- 
particle energies. The dispersion relation ik is pre- 
sented in Fig. |U We assume the unrenormalized single- 
particle energies to follow a tight binding dispersion 
Ek = — 2/j(cosfc — 1), where we set the bottom of the 
band equal to zero and t — 1/2. We recognise two basic 
features for ik- The first is a broadening of the band. 
The effective band width is enlarged compared with the 
original band width W = At. The other one is a decreas- 
ing density of states at the Fermi momentum k = k c F 
and at k = ir — kp. This property is mainly due to the 
dominant peak structure in the spin correlation function 
S(q) at q — 2kp. The wave vector q = 2k p connects the 
two points of the Fermi surface. Therefore the energies 
near the Fermi surface become more strongly renormal- 
ized than energies near the band edge. The pseudo-gap 
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FIG. 4: Renormalized electronic single particle energies it at 
n c = 2/3 and different values J/t. 



like structure at k = kp is thus due to the strong spin 
fluctuations at q = 2k p. As we are going to see this 
behaviour shall have consequences for the electronic den- 
sity of states p(u) . As the Luttinger liquid theory expects 
p(uj) to vanish at oj = a decreasing density of states in 
the renormalized electron spectrum is reasonable. In or- 
der to resolve the observed structures in the renormalized 
electron spectrum we need to examine larger systems. 



B. Static properties 

Let us now study the static correlation functions cal- 
culated in the last section. The first quantity we want 
to consider is the momentum distribution function n(k) 
which is shown in Fig. [5j We obtain meaningful re- 
sults only for couplings up to J/t w 1. This signals a 
breakdown of the flow equation treatment. In order to 
get better results for larger ratios J/t we need to go be- 
yond the third order corrections in the flow equations. 
Looking at n(k) we notice that it is smeared out around 
kp = ti c tt/2. However, we can not decide whether these 
results support the expected Luttinger liquid picture or 
not. The special behaviour of n(k) at k = kp may be 
due to the dominant peak structure of S{q). Since it is 
difficult to resolve the sharp peak of S(q) appropriately 
for a finite system, we are not sure whether the artifact at 
k = kp is due to the finite system size or the approxima- 
tions. Nevertheless, the shape of the momentum distri- 
bution function tends to support a small Fermi surface, 
because there is no feature at k = (n c + \)tt/2. Addi- 
tionally one may question if the effective model (fT2"|) is 
capable of describing a large Fermi surface. The system 
of conduction electrons within the effective model has a 
Fermi momentum n c ir/2. Therefore a singularity in the 
momentum distribution function is likely to be expected 
only at the point k = kp. 

We can get further information from the charge corre- 
lation function G(q). The results are depicted in Fig. [5] 
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via the relation |2£ 



FIG. 5: Moment distribution function n(k) at n c — 2/3 and 
different values J/t. The lines are drawn by ommiting the 
point at k = kp. 



As we expect for small couplings J/t the function C(q) 
takes the form of a noninteracting one-dimensional elec- 



tron gas with a kink at q = 2k\ 



Increasing J/t 



leads to a cusp-like behaviour of C(q) at q = 2k p. In 
addition, the slope at q = drops with growing interac- 
tion stregth J/t. The results displayed in Fig. [5] agree 
ualitatively with the findings from numerical treatments 
in the examined parameter regime (J/t < 1). This 
supports the Luttinger liquid picture of our description. 
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FIG. 6: Charge correlation function C(q) of the electrons for 
different values J/t and band fillings n c . Circles are drawn 
for comparison, taken from [30| at J/t — 1 and n c — 4/5. 

The charge correlation function gives us the possibility 
to derive the parameter K p of the Luttinger liquid theory. 
This parameter is connected to the slope of C(q) at q = 



Kn 



dC(q) 
dq 



(44) 



The outcome is depicted in Fig. [7] as a function of the 
Kondo coupling J/t. As we have already mentioned be- 
fore the slope of C(q) at q = decreases with growing 
coupling strength (up to the allowed value of J/t < 1). 
For vanishing interaction strength K p — > 1 corresponding 
to a noninteracting electron gas. This can be understood 
from the equation for C(q) given in the appendix. Due 
to the flow equation (f53"| for the parameter ~j kq of the 
electron operator transformation all terms vanish which 
represent corrections to the charge correlation function 
of independent electrons. Our findings arc in qualitative 
agreement with recent numerical results from DMRG cal- 
culations [30( | . Xavier and Miranda find a minimum of 
K p (J) at J/t sa 1.5. Remember that our largest possible 
coupling is smaller than 1.5. Quantitatively our results 
are always considerably smaller than the values found in 
ref. [3(| • Another work by Shibata et al. [1] gives results 
for large J/t. In contrast to our findings and to those of 



Xavier and Miranda 
the limit J/t -> 0. 
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FIG. 7: Parameter K p as function of J/t for two different band 
fillings n c . Circles are drawn for comparison, data taken from 
[U for n c = 4/5. 

We also considered the dependence of K p on the band 
filling n c . This is depicted in Fig. [51 The lower possi- 
ble value of the band filling is n c — 1/2 as we do not 
obtain a solution of the flow equations below this value 
within the present approach. At small values J/t we find 
a monotonic decrease by lowering n c . Again we find qual- 
itative agreement with Xavier and Miranda [3O] - As we 
already mentioned in the last discussion our values for 
K p are considerably smaller compared to the numerical 
data. For larger J/t the behaviour deviates even qual- 
itatively from the numerical DMRG data. Whereas in 
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ref. [30] for all values of J/t a monotonic increase was 
obtained on increasing n c , we find a maximum in the 
function K p (n c ). 
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FIG. 8: Parameter K p as function of the band filling n c for 
different values J/t. Circles are drawn for comparison, data 
taken from [30[ for J/t = 0.35. 

The magnetic properties of the one-dimensional KLM 
are significant for the determination of the phase transi- 
tion from the paramagnetic metallic phase into the fer- 
romagnetic phase. The spin correlation function for the 
conduction electrons S cc (q) as well as for the local mo- 
ments Sff(q) show a characteristic increase of the ferro- 
magnetic component q = on approaching the quantum 
phase transition. 
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FIG. 9: Spin correlation function S cc (q) of the electrons for 
different values J/t and band fillings n c . Circles are drawn for 
comparison, data taken from [|| for J/t = 1 and n c = 7/10. 

The spin correlation function of the electrons S cc (q) 
is shown in Fig. [51 The strong peak at q = 2k c F results 



from the sharp maximum of the spin correlation function 
S{q) of the effective model. This can easily be seen from 
the expression of S cc {q) given in the appendix. Another 
characteristic is the finite weight of the ferromagnetic 
component q = which is directly connected with the oc- 
curence of the quantum phase transition. On approach- 
ing the critical J/t the maximum of S cc (q) at q = 2k F 
loses weight in favour of the ferromagnetic component. 
This behaviour marks the phase transition (ij. As we 
already mentioned the present approach is restricted to 
values of J/t < 1. These values are too small compared 
to the value at the transition point which is J/t < 2.5 
for n c — 2/3 [6]. Nevertheless, we observe some tendency 
towards the magnetic phase transition. As in the case of 
the charge correlation function we compare our results 
with numerical data from [8j]. One can clearly see the 
qualitative agreement between the two approaches, al- 
though our findings tend to be smaller than the DMRG 
results. This is important if one considers the points at 
q = 0. The increase of the ferromagnetic component, 
which signals the tendency towards the quantum phase 
transition, turns out to be comparably weak. 
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FIG. 10: Spin correlation function Sjf(q) of the local mo- 
ments at n c = 3/5 and different values J/t 



The situation we have just described is also charac- 
teristic for the spin correlation function Sft(q) of the 
local moments which is drawn in Fig. 1101 For small 
couplings we see that the corrections to the spin corre- 
lation function of the effective model S(q) are negligibly 
small. Even for larger values of J/t we find only small 
corrections. The vicinity of the ferromagnetic component 
q = is shown in the inset. Nevertheless, the qualitative 
behaviour is once again in agreement with numerical re- 
sults [H though the values are somewhat larger. Again, 
the ferromagnetic component gets an increasing weight 
while the q — 2k F component is suppressed. 
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C. Dynamic properties 

In the last section we have presented the results for 
static expectation values and correlation functions. We 
have found that our results are in qualitative agreement 
with numerical data for small couplings J/t. The flow 
equation method sets us in the position to calculate 
not only static but also dynamic correlation functions. 
Within our approach the dynamics of the KLM is de- 
scribed in terms of the effective model flT2j) . Since H is 
blockdiagonal the dynamics for electrons and local spin 
moments seperate. The SBMFT allows us, at least ap- 
proximately, to characterise the excitations of the spin 
system. The excitations of the KLM are determined by 
a noninteracting Fermi gas (conduction electrons) and 
the Schwinger bosons. In this section we shall add new 
aspects to the results obtained by [l(|. 
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FIG. 11: Dynamic correlation functions A+(k,ui) (black) and 
A-(k, uj) (blue) for k = (up) to k — n (down). 

We start with the dynamic properties of the conduc- 
tion electrons. The first quantities we want to consider 
are the electronic spectral functions A±(k,uS) which can 
be measured in XPS and inverse XPS experiments. The 
outcome is shown in Fig. 1111 The energy is measured 
with respect to the Fermi-energy of the conduction elec- 
trons e c F = 7i c 7r/2. As we have already mentioned in 
the last section both functions consist of two parts. A 
coherent quasiparticle-like contribution embodied by the 
finite peak which has a weight Sfc. Its position is sim- 
ply given by the renormalized single-particle energies 
The incoherent background contains pairs of elementary 
excitations of the spin system of H. This follows directly 
from Eq. (|29|) since within the Schwinger boson approach 
for the effective spin system the corrections to the spec- 
tral functions A±(k,ui) are always connected to the cre- 
ation (annihilation) of pairs of bosons. The coupling to 
the continuum of Schwinger boson excitations and the 
results for "fkq give rise to the two maxima around the 
quasiparticle-like peak. 

The electronic density of states p(uj) is an important 
quantity which shows a characteristic behaviour for Lut- 
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FIG. 12: Elektronic density of states p(ui) for J/t — 1 and 
n c = 3/5. 



tinger liquids. It is drawn in Fig. Q21 Again, the en- 
ergy is measured with respect to e c F — n c n/2. We find 
two minima, one in the vicinity of the Fermi level of the 
conduction electrons, ui = 0, the other above the Fermi 
level. This behaviour follows from the pseudogap-like 
behaviour of the renormalized single-particle energies . 
In contrast to our findings DMRG studies from Shibata 
and Tsunetsugu [To| do not yield a minimum but rather 
a peak structure just below uj = indicating the de- 
velopment of a pseudo gap. However, their results were 
performed at finite temperatures. The Luttinger liquid 
theory predicts a density of states following p(uj) ~ \uj\ a , 
< a < 1 in the vicinity of uj — 0. As can be seen from 
Fig. [T2]there is no real vanishing of p(u>) at oj = 0. As we 
are dealing here with a finite system size we are not able 
to resolve p{uj) near uj = and to verify the expected 
behaviour. 

Let us now turn to the magnetic properties. We want 
to present the results for the dynamic spin structure fac- 
tors of the electrons S cc (q,uj) and the local moments 
Sff(q,u>). They describe the magnetic excitations of the 
coupled system and can be measured by inelastic neutron 
scattering experiments. 

We begin with the electronic dynamic spin structure 
factor S cc (q, uj) which consists of a low- and a high-energy 
part. Both are discussed seperately. The low energy sec- 
tor, left panel of Fig. [T3l is characterised by the spin part 
of the effective model TLs, i-e. the continuum of pair ex- 
citations of the Schwinger bosons. The dominant contri- 
bution is therefore found at q = 2k c F . It is multiplied by a 
factor 1/4 for a better comparison. We also see that there 
are regions where no excitations are possible. Further- 
more, the gap in the spectrum of the elementary excita- 
tions Cj q leads to a gap in the low-energy part of S cc (q, uj). 
The high-energy sector of S cc (q, uj) is shown in the right 
panel of Fig. [13] The spectral weights are about 10 times 
smaller compared to the weights of the low-energy part. 
The main contribution arises from electronic particle-hole 
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excitations of the effective model. The specific form of 
this contribution shows therefore the characteristics of 
a one-dimensional electron gas: a gapless excitation at 
q = 2kp and regions between < q < 2k F where no ex- 
citations are possible. In addition to the terms describing 
pure particle-hole excitations there are also terms involv- 
ing the elementary excitations Cu q of the effective spin 
system. They are responsible for the broadening of the 
structures in the high energy sector of S cc (q,u>). 
The DMRG calculations of Shibata et al. for S cc (ui) 
j ^S cc (q,uj) showed a small peak at very low energies 
and a larger double peak structure at higher energies [l(| • 
We obtain a similar peak structure, but in contrast to the 
results of [10J the spectral weight of the low energy part 
is much larger than the spectral weight of the high energy 
part. This does not agree with the picture of an exhaus- 
tion of the electronic low-energy spin degrees of freedom 
due to singlet formation described by [10| ■ 

Finally we want to discuss the magnetic excitations 
of the system of local spin moments described by the 
dynamic spin structure factor Sff(q, ui). As in the case of 
the electronic spin structure factor S cc (q, ui) this function 
comprises a low- and a high-energy part. The first one 
is again determined by the elementary excitations Cj q of 
the spin part of the effective model Hs- It is depicted in 
Fig. Q3] and possesses the same features as the low-energy 
part of S cc (q, ui). From this picture we can clearly see the 
influence of the low-energy spin excitations. The distinct 
structure at q = 2kp gives rise to the pronounced peak 
in the static spin correlation function S(q). Once again 
we point out that the energy scale of these excitations 
are quite small compared with the effective band width 
of the electrons. The DMRG results of Shibata et al. for 
Sff(oj) = J ^Sffjq, lo) show a large peak structure at 
very small energies [10]. They assume that this is due 
to collective spin excitations of the Luttinger-liquid. In 
our approach the low-energy peak is the result of the 
continuum of elementary excitations of the effective spin 
system, which we described in terms of Schwinger bosons. 



FIG. 14: Low energy part of the dynamic spin structure factor 
Sff(q,u>) of the local moments for J/t = 1.5 and n c = 3/5. 
The colour function is normalised to the maximum contribu- 
tion. 
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FIG. 15: High energy part of the dynamic spin structure fac- 
tor Sff(q,ui) of the local moments at J/t = 1.5 and n c = 3/5 
for q = (bottom) and q — it (top). 

The high-energy part of Sff(q, ui) is shown in Fig. [T5] 
As in the case of S cc (q,u>) it is mainly determined by 
particle-hole excitations of the Fermi sea. At larger cou- 
plings J/t the elementary spin excitations Co q lead to the 
broadening of the peak structure. Shibata et al. obtain a 
second peak in the high-energy sector of S//(w) [lQJ. Our 
approach yields a similar structure in the local spin dy- 
namics, although the spectral weight of the high-energy 
part is much smaller than the spectral weight of the low- 
energy part. 

We further note that Sf/(q, ui) exhibits a finite gap. This 
is an artifact and due to the approximations that we have 
made for the spin operator transformation. By taking 
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into account higher correction terms the transformed spin 
operator of the local spin moment couples to the spin op- 
erator of the conduction electrons. This gives rise to a 
gapless mode in Sff(q,u>). 

V. CONCLUSION 

In summary, we used the method of continuous uni- 
tary transformations (flow equation method) to examine 
the one-dimensional KLM. The renormalization proce- 
dure was employed to integrate out the coupling between 
conduction electrons and local spin operators. In that 
way we derived an effective Hamiltonian which consists 
of an one-dimensional noninteracting electron gas and a 
Heisenberg chain interacting via an RKKY-like coupling. 
In order to treat the spin chain we used a Schwinger bo- 
son mean field theory (SBMFT). Thereby we were able to 
calculate static and dynamic correlation functions. The 
investigation of the electronic momentum distribution re- 
vealed a small Fermi surface. We gave arguments, re- 
ferring to the effective model, why we were not able to 
obtain the large Fermi surface scenario in our approach. 
Nevertheless, the static spin and charge correlation func- 
tions of the electrons agreed qualitatively with numerical 
results. In addition we obtained the parameter K p of 
the Luttinger liquid theory and found also qualitative 
agreement with recent DMRG calculations. The present 
approach was restricted to parameter regimes J/t < 1- 
Although the quantum phase transition from the param- 
agnetic metallic into the ferromagnetic phase takes place 
at larger values, we observed some tendency to a stronger 
ferromagnetic component in the static spin correlation 
functions. The new aspect of this work was the extension 
of calculations for dynamic properties by means of the 
flow equation's method. We showed that the electronic 
spectral functions comprised a coherent quasiparticle-like 



peak determined by the renormalized electronic disper- 
sion relation. The coupling to the low-energy excita- 
tion of the effective spin model gave an incoherent back- 
ground comprising two maxima near the quasiparticle- 
like peaks. Finally, we also computed the magnetic ex- 
citations of both the electrons and the local spins. The 
corresponding spin structure factors always consisted of a 
low-energy part, determined by the Schwinger boson pair 
excitations, and a high-energy part, mostly determined 
by electronic particle-hole excitations. The latter there- 
fore showed the special features of the one-dimensional 
Fermi surface. The electronic spin structure factor exhib- 
ited a gapless mode at q = 2k F . Our results for the elec- 
tronic spin dynamics did not agree with the exhaustion 
picture described by [1(3]. The gapless mode at q — 2k p 
should also be seen in the spin structure factor of the lo- 
cal moments. There we argued that further corrections in 
the spin operator transformation would lead to a gapless 
mode. 
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APPENDIX A: CORRELATION FUNCTIONS IN 
THE FLOW EQUATION APPROACH 

In this appendix we give the rather lengthy expressions 
for the static and dynamic correlation functions omitted 
in the text. These are the static charge correlation func- 
tion C(q) which reads 



N 

kk' aa' 

2 

k 

+ Jp Y 5fe 5fe +9 >< fc +<? %-i> k S ( k -P-l) n k+ q {^ ~ n k ) + 7k-q, p S(k -p-q)n k (l- n p ) 

kp kp 

+ J^pY & klk+ q ,p S(k - p + q) n p (l - n k ) + & k&p%+q,k%-g, P S(k-p- q) n k {\ - n p ) (Al) 

kp kp 



The dynamic spin structure factor S cc (q,cj) of the electrons takes the form 

/oo 
dt (s,(t) • s_ 9 > e Mt 
-oo 



2 ^ 



2N 
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CKfe % lp+q,k %-q,p ( u p' v p'+k-p-q ~ Vp'Up' +k-p-qf n k(^ ~ n p )8(uJ - Q p t - UJ p > +k-p-q ~ ?k + £p) 



kpp' 

4N 3 



+ AM3 ^2 lk+q,p (Up'Vp'+k-p+q ~ V U P' + k-p+q) 2 M 1 - n k )S(u - Ul p , - U)p' + k-p+q ~ £p + £k) 
kpp' 

J^3 ^2 ^ 7k-q,p ( U p' V p' + k-p-q - Vp'Up'+k-p-q) 2 U k (l - n p )5(uj - U) p > - LUp' + k-p-q - £k + Ep) 
kpp' 



+ 2^y3 ^2 Oik a p (jk+q,k %-q,p + Jk~q,k lp+q,p)nk U p (Up,V p , +q - Vp,Up, +q ) 2 S(UJ - U) p > - Cjp'+q) 
kpp' 

+ -jjy^ ^2 ® k ®P (7*1-9,* 7p-g,p + lk+q,k %+q,p)nk Up (u p 'V p/+q ~ Vp'Up' +q ) 2 S(U! - U! p < - UJp'+q) (A2) 
kpp' 

Here, it can clearly be seen that the second and third line involves only particle hole excitations of the Fermi sea of the 
effective model. The last two lines represent the low energy sector of S cc (q, u>) as they include only pair excitations of 
Schwinger bosons. On integrating over the energy to one obtains the expression for the static spin correlation function 
Scc(q). ^ 
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